Computing and Sampling Restricted Vertex 
Degree Subgraphs and Hamiltonian Cycles 

by Scott Sheffield 

Abstract 

Let G = (V, E) be a bipartite graph embedded in a plane (or n-holed 
torus). Two subgraphs of G differ by a Z -transformation if their symmet- 
ric difference consists of the boundary edges of a single face — and if each 
subgraph contains an alternating set of the edges of that face. For a given 
(j) \ V i — * Z + , S<f> is the set of subgraphs of G in which each v 6 V has 
degree (j>(v). Two elements of are said to be adjacent if they differ by 
a ^-transformation. We determine the connected components of S$ and 
assign a height function to each of its elements. 

If <f) is identically two, and G is a grid graph, S$ contains the partitions 
of the vertices of G into cycles. We prove that we can always apply a series 
of ^-transformations to decrease the total number of cycles provided there 
is enough "slack" in the corresponding height function. This allows us to 
determine in polynomial time the minimal number of cycles into which G 
can be partitioned provided G has a limited number of non-square faces. 
In particular, we determine the Hamiltonicity of polyomino graphs in 
0(|V| 2 ) steps. The algorithm extends to n-holed-torus-embedded graphs 
that have grid-like properties. We also provide Markov chains for sampling 
and approximately counting the Hamiltonian cycles of G. 

1 Introduction 

Let G — (V,E) be a planar or n-holed-torus-embedded polygon graph (i.e., a 
graph whose faces are simple polygons) and let 6 be a function on the oriented 
edges of G. Then the height function space % = (G, 8) is the set of (real or 
integer valued) functions 4> on the face set F of G satisfying difference restrictions 
4>{f2) — <A(/i) < S(e) whenever e is an edge shared by fa and fa oriented with 
fa on the left. 

If G is planar, "K is in one-to-one correspondence with the set of closed elec- 
tric flows u> on G satisfying — 6(e~ 1 < w(e) < 6(e) for each oriented edge e. This 
allows us to convert questions about randomly chosen network flows or corre- 
sponding objects (e.g., perfect matchings on bipartite graphs) into questions 
about random height functions. 

One advantage of this conversion is that recently introduced perfect sampling 
techniques, including Coupling from the Past (see Q) and Fill's algorithm (see 
[Ejl), apply especially well to height function spaces. Another is that large 



deviations principles and asymptotic limiting shapes can sometimes be derived 
for random height function spaces in a way that sheds light on local behavior. 
(See § and @ for numerous references to examples involving domino tilings.) 

This paper uses height function theory to solve two other problems that 
have been dealt with in the literature. The first of these is Z -transformation 
connectedness. The problem here is as follows. Let G be a bipartite, planar (or 
n-holed torus embeddable) graph G = (V,E), and <fi a function <fi : V i— > Z + . 
Let S,/, be the set of all subgraphs H of G satisfying degn{v) — 4>{v) at every 
vertex v. 

Two subgraphs of G are said to differ by a Z-transformation if their sym- 
metric difference consists of the boundary edges of a single cell — and if each 
subgraph contains an alternating set of the boundary edges of that cell. We say 
two elements of are adjacent if they differ by a single Z-transformation. The 
problem is to understand the connected components of S^. 

In the case that <f> is identically one (and thus is the set of perfect match- 
ings of G) the problems has been solved explicity for polyomino graphs, simply 
connected hexagon graphs, and other families without the aid of height functions 
(see II, ||, @, @ and g|). 

Using height functions, however, we can always determine the structure of 
S,/, and enumerate its connected components. The key is to observe that two 
subgraphs differ by a Z-transformation if and only if their corresponding height 
functions are equal everywhere except on a single face where they differ by 
one. In fact, we will extend our definition to say that any two height functions 
differ by a Z-transformation if this is the case. Given this formulation, the 
connectivity results are almost trivial. 

Similarly, several papers (including those in @, 0, 0, g], and @) 
have used a variety of techniques to produce Hamiltonian cycles in grid graphs 
without using height functions. However, these algorithms and characterizations 
do not apply in general. For example, ]22| solves the Hamiltonian path problem 
for rectangular grids; 0] and JjT] deal with rectangular grids with one or two 
vertices removed; {h| defines a special class of Hamiltonian grid graphs by 
inductively adjoining pairs of rows; and || gives a linear algorithm for the 
Hamiltonicity of grid graphs shaped like staircases (with possibly irregular step 
sizes). 

Using height functions, we can go much further. First, we specialize to the 
case that G is a grid graph and (j) is identically two (so that elements of S$ are 
cycle covers of G). We will see that applying a series of Z-transformations to 
lattice square faces sometimes allows us to join those cycles together. This turns 
out to always be possible provided there is enough "slack" in the corresponding 
height functions; our algorithm will search through the connected components 
of 8$ until it finds one for which the corresponding height functions do have 
enough "slack." 

We will always assume that G is 2-connected (otherwise, it is clearly not 
Hamiltonian), and thus, the faces of G are simple polygons without repeated 
edges. The complexity of our algorithms will depend on m, the number of faces 
of G that are not lattice squares. (Intuitively, these correspond to "holes" in the 
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grid graph.) We give a (probably not sharp) upper bound on the runtime by 
0(|F| 2+ "2"), where F is the set of square faces of G. In particular our algorithm is 
Od-Fl 2 ) for grid graphs with no non-square faces (also called polyomino graphs). 
Similar results hold for Hamiltonian paths with fixed endpoints and for other 
gridlike families of graphs embedded in n-holed tori. 

An algorithm very similar to the one we present was published in ]39| ]; this 
paper (which applies only to polyomino graphs and gives a somewhat weaker 
bound of 0(|y| 4 )) does not discuss height functions explicitly, but it contains 
many definitions and constructions that are most naturally understood in height 
function language. For example, the somewhat unintuitive definition of the dis- 
tance between 2-factors in ]39| (based on an induction using nestings of "alter- 
nating cycles") is the l\ distance between the corresponding height functions. 

An intuitive reason we might expect the Hamiltonicity problem to be more 
difficult for grids with holes is that the path could belong to any one of the ex- 
ponentially many homotopy classes of non-intersecting cycles in the plane minus 
the holes. An interesting open problem we don't solve is whether, given a homo- 
topy class, one can determine in polynomial time whether there is a Hamiltonian 
cycle in that homotopy class. To compensate, we provide a related result that 
is neither strictly stronger nor strictly weaker: we produce a correspondence 
between cycle partitions and network flows and provide a polynomial algorithm 
for the existence of a Hamiltonian cycle whose corresponding network flow lies 
in a particular homology class of the plane minus the holes. 

Our last topic will be the enumeration, sampling, and counting of the Hamil- 
tonian cycles (or Hamiltonian paths) for a given grid graph using Markov chains. 
Although the convergence rates of the Markov chains are not known in general, 
we can provide heuristics and upper bounds in special cases. 

Finally, as a potential application, we point out that theoretical chemists and 
statistical physicists make frequent use of Hamiltonian paths and cycles in grid 
graphs to model proteins and other heteropolymers. These results may prove 
useful in analyzing the phase transitions and other properties of two-dimensional 
single-chain heteropolymer models. 

2 Height Functions: Definitions and Basic Re- 
sults 

The results in this section are standard material from combinatorial topology 
and the cycle decomposition of planar electric flows; we include them to fix 
notation and make the paper self-contained. 

2.1 Defining Chains and Homology 

From here on, we assume that G is embedded in a plane or n-holed torus and 
we will refer to the connected components of the surface of embedding minus 
G as faces. One easily checks that if G is connected and embedded in a surface 
with genus as small as possible, all of the faces are simply connected and have a 
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sequence of (possibly repeated) oriented edges e\, . . . , as a clockwise oriented 
boundary. Because of the orientability of our surface, if a face contains the 
same edge twice on its boundary, it must be oriented in opposite directions each 
time. If none of the faces of G has repeated edges on its boundary, then all 
face boundaries are simple polygons, and we say G is a polygon graph. We have 
already defined ^-transformations for polygon graphs, but we can extend that 
definition. If a face / has repeated edges, we say two subgraphs Hi and H2 
differ by a Z-transformation at / if their symmetric difference is the set of non- 
repeated edges of / and if one subgraph contains all non-repeated edges oriented 
from black to white on the clockwise boundary of / and the other contains all 
such edges oriented from white to black. 

Let So(G), Si(G), and 82(G) be the formal vector spaces generated respec- 
tively by the vertices, oriented edges, and faces of G (where in S\(G), an edge 
oriented one direction is the additive inverse of the same edge oriented the op- 
posite direction). For example, if G has five faces, 82(G) is the five-dimensional 
vector space of sums of the form X)/gf a //i with a/ G R. 

The dual perspective is to think of an element of 5*2 (G) as a function from 
these five faces to the real numbers. Because we will need this dual perspective 
frequently, by slight abuse of notation, whenever r is an 2-chain, we will write 
r(/) to mean the value of cif in the sum r = X)/g_f a ffi (i- e -' the inner product 
of t(/) and /). We follow similar conventions for S\(G) and Sq(G). Elements of 
Si(G) are called i- chains. There is an obvious correspondence between 1-chains 
and electric flows. 

Next, we define a boundary operator: 

d : 5i(G) ^ S Q (G) 
d : S 2 (G) ^ S^G) 

as follows; if e — (v, w) is an edge (oriented from vertex v to vertex w), then 
d(e) = w — v. If / is a face, then d(f) is the sum of the edges on the boundary 
of /, clockwise oriented. ("Clockwise" is well-defined because the plane and n- 
holed torus are orientable.) Note that if the boundary of / contains a repeated 
edge, both that edge and its inverse will be included in the boundary sum and 
will cancel each other out. 

If to is a 1-chain, we say ui is closed if d(uj) = 0. We say to is exact if u> = d(r) 
for some chain r. Since, d 2 (f) is clearly equal to zero for any single face /, 
d 2 (r) — for any 2-chain r. Hence, every exact one-chain is closed. 

The homology H\ (G) is defined to be the vector space of closed 1-chains mod 
the space of exact 1-chains. In other words, homology represents the failure of 
the sequence 

1 — > 82(G) 1 — >^ Si(G) 1 — >^ Sq(G) 1 — > 

to be exact. The most important result from combinatorial topology we use is 
the following: 
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Theorem 1 The dimension of the homology space depends only on the topolog- 
ical structure ofT, the surface of embedding, not on the graph G. The homology 
is trivial if T is the plane and 2n- dimensional if T is the n-holed torus. 

The reader who wishes to verify these facts can do so with the following 
steps. 

1. Verify that the dimension of dS 2 (G) is \F\ — 1 by checking that the kernel 
in the space S 2 of the map d is one-dimensional, spanned by J2feF /• 

2. Verify that the dimension of dS\(G) is |V| — 1 by showing that for any 
two vertices, v, w, the 0-chain v — w is in dS\(G) and then showing that 
these elements span a codimension one subspace of S*o(G). 

3. Conclude that the kernel of d in Si(G) has dimension \E\ — \V\ + 1 and 
the homology has dimension 2 — |F| + |S| — \ V\. Compute this value using 
the Euler characteristic. 



2.2 Defining Z-transfor mat ions and in Terms of Chains 

Wc refer to vertices of G as black and white depending on which of the two 
partite classes they belong to. Given a subgraph H of G, we then define the 
1-chain u>h to be the sum of all the edges in H, oriented from black to white. 
Then H e Sj, if and only if 

dujH = ^ (p(v)e(v)v 
veG 

where 




1 v is a black vertex 
— 1 v is a white vertex 



From this it follows also that if Hi and H 2 are both in S4, dujH 1 — dujH 2 — 
0. In other words, — luh 2 1S always closed (though not necessarily exact, 
unless G is planar). Furthermore, two subgraphs Hi and H 2 differ by a Z- 
transformation if and only if — ^h 2 is equal to df or — df for some single 
face /. 

This implies that if one can move from Hi to H 2 by a sequence of Z- 
transformations (i.e., Hi and H 2 are in the same connected component of 
then ojh 1 — ll>h 2 can be written as a sum J2f eF dfdf where a/ are integers. In 
particular, ojh 1 — loh 2 is exact. 

Whenever loh 1 — ujh 2 is exact, we will say that Hi and H 2 are homologically 
equivalent and write Hi ~ H 2 ; since this is clearly an equivalence relation, it 
partitions into homology classes. 
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2.3 Chains and General Network Flow Problems 



In the above formulation, a 1-chain u is equal to ujh for some H e if and 
only if 

1. gLj = 4>{ v ) e ( v ) v 

2. For each black to white oriented edge e, ui(e) 6 {0, 1}. 

Equivalcntly, we might require that u> represent an electrical network flow 
such that each black vertex is a source of <j)(v) units of current, each white vertex 
is a sink of 4>(v) units, and each edge conducts one or zero units from black to 
white. This is best understood as a special case of the following well-known 
constrained network flow problem: 

Given a graph G, an upper bound 6(e) and a lower bound (given by — 5(e^ 1 )) 
on the amount of current passing through each directed edge e, and a function </> 
on the vertices dictating the net amount of current (positive or negative) flowing 
into that vertex, describe the set of (real or integer-valued) flows with this 
property, i.e., the set of one-chains w such that du = ^ <p(v)v and u)(e) < 5(e) 
for each directed edge e. 

The first step is to reduce to the case that <j) is identically zero (and hence 
our one-chains are closed) as follows. Let a be any (integer or real valued) flow 
such that the amount of flow into v is <fr(v) for each v. Instead of seeking u>, 
we can look for ZJ = uo — a such that dZJ = and ZJ(e) < 5(e) for all e, where 
5 = 5 — a. 

Next, we treat the set of closed ZJ one homology class at a time. (If G is 
planar, there is only one homology class.) Accordingly, we restrict ourselves to 
u> such that ZJ is exact. (Different choices of a correspond to different homology 
classes.) Then for each such ZJ, there is a 2-chain t such that dr = ZJ. Viewed in 
the dual sense as a function on the faces of G, t is the height function discussed 
in the introduction, and it is uniquely determined up to a constant. To eliminate 
ambiguity, we will require that r(/o) = for some reference face fa. If f\ and ji 
share an edge e (oriented with /i on the left), the condition that ZJ(e) < 5(e) is 
equivalent to the condition that r(/2) — t(/i) < S( e )- We have now reduced the 
characterization of flows in a homology class to the characterization of height 
functions satisfying difference restrictions. 

To help us better understand the restrictions, for adjacent /i and fi, we 
define d(fi, $2) = 5(e), the largest allowed value of r(/2) — Next, suppose 

/1 and f 2 do not necessarily share an edge, but P = (fi,f a (i),- ■ ■ , f<T( m ),h) is a 
path connecting fi and f 2 (so that each f a u) shares an edge with the face before 
and after it). Then we define D P (f 1 , f 2 ) = d(f 1 ,f a{l) ) + d(f a{l) ,f lj{2) ) + ...+ 
d{fa{m),h)- We then deduce the useful bounds, -D P (f 2 ,f 1 ) < T H (f 2 )-r(f 1 ) < 
Dp(fi, f 2 ). When it exists, we define D(fi, f 2 ) to be the minimum value of 
Dp(fi,f 2 ) as P ranges over all possible paths connecting fx and f 2 . This 
includes single-element paths; thus, if there are no negative cycles (i.e., no paths 
P from a face / to itself such that D P (f, f) < 0), then D(f, f) = 0. 
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Theorem 2 There exists a function r on the faces of G satisfying the difference 
restrictions if and only if there exists no path P from a face f to itself such that 
D P {f,f)<{). 

The above bounds clearly imply that the latter condition is necessary. To 
prove that it is sufficient, we first note that if there are no negative cycles, then 
for any f\ and fi, a minimum Dp(fi, fi) must exist. To see this, note that if 
P' is a path connecting f 2 back to f± and P + P' is the concatenation of P and 
P', then 

< Dp +P ,(fx, h) = D P (f u f 2 ) + D P ,(f 2 , fi) 

If we can make Dp(f\, fi) arbitrarily small, we obtain a contradiction by choos- 
ing P so that D P+P i(f 1 , fi) is negative. 

Now, to prove existence of a r satisfying our restrictions, we show that there 
is one choice of r that achieves its maximal allowed value on every face: that 
is, t(/) = D(fo, f) for all /. To see that this function satisfies our restrictions, 
suppose fi and fi are joined by an edge e with fi on the left. Then one easily 
checks that 

D(fo,h)<D(f ,f 1 )+d(f 1 ,f 2 ) = D(fo,f 1 )+5(e) 
D(f , fi) < D(f , f 2 ) + d{f 2 ,h) = D(f , h) - 1(e) 

and hence 

r(f2) - r(h) < 6(e) 
From this construction of r, we get the following corollary. 

Corollary 1 If a reference face fo is required to have t(/o) = 0, then there 
exists a r satisfying our difference restrictions such that t achieves its maximal 
(or symmetrically, minimal) possible value on each face of G. In other words, 
if we write T\ > t% whenever T\(f) > T2(/) for all f, the resulting poset of 
allowable height functions, contains maximal and minimal elements. 

The distance function D can be thought of as representing a (non-symmetric 
and possibly negative) cost of traveling from f\ to fi and can be efficiently 
computed by standard shortest path and negative cycle detection algorithms, 
such as Dijkstra's algorithm (see 0) and its many variants. 

Clearly, every homology class of S\(G) corresponds to the class of height 
functions arising from some choice of a. Now, fix ceo to be any chain satisfying 
dcto = 4>( v ) v > an d let ci, . . . , c„ be a generating basis of the elements of the 
homology space of G; (such a basis can be computed with simple linear algebra, 
or Z- module reduction in the integer- valued case). We can now enumerate 
the homology classes by determining for what values of a = (ai, . . . ,a n ) the 
one-chain a = a + a^Cj gives rise to negative cycles in the corresponding 
restricted difference problem. 

It is clearly enough to consider non-self-intersecting paths P, of which there 
are finitely many. Given a path P from a face / to itself, note that the value 
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Dp(f, f) has affine dependence on a, and thus the restriction that Dp(f,f) > 
forces a to lie one side of a hypcrplanc. Hence, for the real (integer-valued) flow 
problem the allowable o's correspond to real (integer-valued) vectors inside a 
convex polytope. Furthermore, for any allowable a and a' and any i, one easily 
checks that \cn — a\\ is bounded above by |-F|M, where M is the maximum size 
of a restricted difference interval. Thus, in the integer-valued case, (|.F|M)™ 
is an upper bound on the number of homology classes. In particular, for the 
restricted vertex degree problem, \F\ n is a (usually very weak) upper bound for 
the number of homology classes of S$. 

3 Connectedness Results 

3.1 Connected Components of Height Function Spaces 

In this section, we study the connectedness of the set of flows in a given homology 
class; recall that once the homology class is given, allowable flows are in one to 
one correspondence with functions r on the faces of G such that whenever /i 
and f 2 share an edge e (oriented with /i on the left), 7(e) < t(/ 2 ) — r(/i) < u(e). 

Now, given a set of faces and difference restrictions, we say two faces /i 
and f 2 are in the same face cluster (or f\ ~ f 2 ) if D(fi,f 2 ) = —D(f 2 ,fi). 
Equivalently, f\ ~ / 2 if for some path P from f\ to itself which passes through 
/ 2 , Dp(/i, /i) = 0. This is clearly an equivalence relationship, and moreover, if 
/i ~ f 2 , one easily checks that D(f 1 ,f 3 ) = D(f u f 2 ) + D(f 2 , / 3 ) and D(/ 3 , /i) = 
D(fa,f 2 ) + D(f 2 ,fi) for any third face / 3 . Whenever fx ~ f 2 , the difference 
r(/i) — r(f 2 ) is fixed to be exactly D(f\, / 2 ) for any allowable height function; 
thus, it is always impossible to apply ^-transformations to an allowable r at fx 
or f 2 . 

The following result is now easy to prove, and is based on the intuitive 
notion from linear programming that if a convex polytope is full-dimensional 
in its space of embedding, any two points (or at least interior points) can be 
connected by paths that only move in one coordinate direction at a time. 

Theorem 3 The Z -transformation graph on 5^ is connected if and only if there 
is at most one face cluster containing more than one face. If there are k face 
clusters (k > 1) with representatives /o, ■ • • ,fk-i> then (taking f as the ref- 
erence face so t(/o) = 0), the connected components of the Z -transformation 
graph are in one to one correspondence, with the set of values t(/i), . . . , r(fk-i) 
satisfying -D{f j ,f i ) < r{fj) - r(fi) < D(fi, fj) for all < i, j, < k - 1. 

First, since we cannot apply ^-transformations to make changes to t(/) 
when / is on a fixed face, any n and t 2 are clearly in different connected 
components if they disagree on any of the fi, ■ ■ ■ , fk-i- 

For the other direction, we will need some definitions. For each path P from 
a face / to itself that includes at least one face other than /, define the length \P\ 
to be the number of edges of the path (so (/, fi, f 2 , f) has length three). Then let 
(3 be the minimal value of Dp(f, f)/\P\ as / ranges over all faces and P ranges 
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over all non-trivial paths from / to / that do not intersect themselves except 
at their endpoints. Next, define the distance \t\ — t 2 | = X)/gf l r i(/) — T 2(/)|- 
We claim that with at most Tl ~ r2 + \F\ Z-transformations, we can transform 
T\ into r 2 . To see this, first let 7 be the maximal value assumed by Ti(/) — t 2 (/) 
as / ranges over the faces in F. We can assume without loss of generality that 
this value is positive, and / is a face on which this value is achieved. Next, let 
m be the minimal value such that decreasing the value of n(/) produces an 
allowable height function. Then if m > t\ (/) — t 2 (/) , we decrease the value of 
T i(.f) by ri(/) — i~2(f), making the two height functions equal at that face. If 
this is not the case but m > /?, we decrease the value of t\ (/) by m, reducing the 
distance |ti — t 2 | by at least (3. Finally, suppose m < ti(/) — r 2 (/) and m < (3. 
Then there must be a neighboring f\ such that decreasing Ti(/i) by m causes 
r 2 — n to be greater than ci(/i,/ 2 ). In other words, there is a neighboring fx 
such that, d(f, /1) — (ti(/i) — ri(/)) < m < [3. Next, we go through the above 
steps again, either reducing the distance by at least (3, creating equality on one 
face, or finding an / 2 with d(fi, / 2 ) — (ti(/2) — T i{fi)) < m < (3. The sequence 
/> /i) /2j ■ ■ ■ we produce here can never intersect itself, because if . . . , /; were 
are cycle P formed by this process, we would have 

1-1 1-1 
D P (f, /) = E d (/*> - ( T i(^) - - E /i+i) < (* - 

and this contradicts our assumption about the minimality of (3. Hence, we can 
cither decrease the distance by (3 or produce equality on a face at every step, 
and we can apply a similar procedure for the faces on which n(/) < T 2 (.f); 
from this we deduce the bound on the number of steps. Finally, to see that 
there exists at least one possible height function for each allowable set of values 
r(/i), . . . , r(fk-i), we can use an extrapolation argument to set each r(/) equal 
to its maximal possible value. That is, r(/) is the maximum of r(fi) + D{fi, /) 
as i ranges from to k — 1. The proof that this yields a valid height function is 
analogous to that given earlier for height function existence. □ 

In the integer-valued case, every such step decreases the distance between 
Ti and t 2 by at least one, so we get the following corollary. 

Corollary 2 In the case of integer-valued network flows, if T\ and t 2 agree on 
all of the face clusters with more than one element, T\ and r 2 can be connected 
by at most t 2 — t\\ Z -transformations. 

Furthermore, when we are interested in random sampling or Coupling from 
the Past, we can always make our space of flows connected by allowing Z- 
transformations on face clusters instead of just individual faces. Further modi- 
fications may be made to increase the coupling speed in many cases, but these 
are beyond the scope of this paper. 
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3.2 Restricted Vertex Degree Subgraph Connectedness 

In the Restricted Vertex Degree Subgraph case, we can restate the connectedness 
result using the following language. Wc say an edge of G is free if there exist 
Hi,H 2 £ Sfj, such that the e is contained in Hi and not H 2 . Otherwise, we say 
the edge is fixed; i.e., it is necessarily included (not included) in each element 
of Sq. Similarly, we say an edge e of G is homology class free with respect to a 
given homology class if there exist Hi , H 2 in that homology class such that e is 
contained in Hi and not H 2 ; otherwise, e is homology class fixed. Clearly, two 
adjacent faces are in the same face cluster if and only if the edge between them 
is homology class fixed-and in the planar case, homology class fixed edges and 
fixed edges are equivalent. 

The notion of free and fixed edges is used in organic chemistry. If we treat 
carbon atoms as vertices, single or double bonds as edges in G, and double bonds 
as edges in H, benzene molecules exist only when G is a hexagon graph that 
admits a perfect matching; when multiple matchings are present, the double 
bond configuration occupies a quantum superposition of possible matchings, 
and carbon atoms joined by free edges of G have different bond lengths from 
those joined by fixed edges. 

Using this definition, we have another formulation of our connectedness re- 
sults: 

Theorem 4 The following are equivalent: 

1. There exists a 1-chain a defined on G such that da = ^2<j)(v)e(v)v and 
each edge e oriented from black to white satisfies < a(e) < 1. 

2. There are no fixed edges in G. 

To sec that the first implies the second, first, as in previous sections, we 
use a to define ujh = Wi? — a and we choose th to be the 2-chain such that 
Th(/o) = for a reference face /o and drn = u>h- Note that the statement 
implies that for any pair of adjacent faces, fi and f 2 , < |d(/i,/2)| < 1- In 
particular, this implies that for any non-trivial path P from a face / to itself, 
Dp(f, f) is positive, and hence each face cluster contains only one face. 

For the converse, if there are no fixed edges on G, then let Hi,.. . ,H^ be 
some set of subgraphs in S$ such that edge is contained and not contained in 
at least one subgraph. Then a — i satisfies the requirements. □ 

Note that if we let if i, ... , Hk run through all the subgraphs in Sq, then 
a £ and hence j ^ = 0. We then conclude that £ T ff; = 0; because 

it has zero as its boundary and because t^(/o) is always equal to zero. This 
is the unique choice of a such that the average of the corresponding height 
functions is zero. Clearly, the problem of determining this a — the average 
value of uiHi — is equivalent to the problem of determining the average height 
function when a different, fixed a is chosen. 

Using the same proof, we can generalize to clusters of faces. 

Theorem 5 The following are equivalent: 
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1. There exists a 1-chain a defined on G such that da = ^ 4>(v)e(v)v and for 
a partition of the set of faces F\, . . . ,Fk, each edge e separating members 
of different Fi oriented from black to white satisfies < a(e) < 1. 

2. No edge separating members of different Fi is fixed. 

This follows from the fact that if an edge e were fixed, then we would have 
fi ~ fi for the two opposing faces and hence there would exist a closed path P 
starting and ending at f\ and passing through fi with Dp(/i, f\) = 0. However, 
we know that all the terms in the sum for Dp(f\, f\) are equal to or greater 
than zero and at least one (one separating an element of the first Fi from the 
second) is positive. □ 

Similarly, we can state the result for particular homology classes. This latter 
formulation is useful because the second condition implies that all subgraphs of 
G corresponding to the homology class defined by a (i.e., all H such that loh — ao 
is exact) can be connected by a series of Z-transformations. 

Theorem 6 The following are equivalent: 

1. There exists a 1-chain a defined on G such that da = ^2 4>{ v ) e { v ) v > a — a a 
is exact, and for a partition of the set of faces Fi,... ,F^, each edge e 
separating members of different Fi oriented from black to white satisfies 
< a(e) < 1. 

2. No edge outside of £ is homology class fixed with respect to the homology 
class of H for which u>h — ao is exact. In other words, no faces from 
different Fi are part of the same face cluster. 

Finally, if fo is the unbounded face of a planar graph G, we will sometimes 
be given a natural choice of a such that a — ll>h is closed at all interior vertices 
of G but not necessarily at those vertices V\, . . . which border on / ; in this 
case, our height functions cannot be well defined on / . However, suppose we 
add a new vertex v to the middle of f together with an edge from v to each 
vertex on the boundary of fo, thus dividing fo into triangles, one triangle for 
each edge. Then for each edge connecting v to Vi, we define a(e.i) so that 
a — ujH a is closed at Vi. Having done this, it follows that a is also closed at v 
(since v cannot be a source or sink if no other vertex is a source or a sink) , and 
our height functions th are now well-defined for any H in S^; and if we choose 
one of the triangles in the division of / as a reference face, we see that the value 
of th on each of these triangles is independent of the choice of H . We will refer 
to these as the boundary values, noting that one such value corresponds to each 
edge of Jo- 

One reason for adopting the above perspective is to force a to satisfy some 
regularity conditions. For example, suppose G is a bipartite planar graph such 
that every vertex not on that unbounded face has exactly has k edges and that 
<j> is identically equal to j (with < j < k). Then can we set a identically 
equal to j/k on each edge oriented from black to white. For example, in the 
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case of perfect matchings on a polyomino graph, we can take a(e) to be 1/4 
everywhere for each edge oriented from black to white. This is essentially the 
approach used to describe height functions in, for example, M. Similarly, if G 
is a hexagonal graph with no "holes" (i.e., all unbounded faces are hexagons), 
we can understand the perfect matchings of G by taking a — 1/3. Applying 
the connectedness arguments of the previous section, we can also deduce the 
following: 

Theorem 7 Let G be a bipartite planar graph such that every vertex not on the 
unbounded face has exactly k edges and 4> is constant. Then is connected. 

In the case of partitions of polyomino graphs into cycles, we can take a(e) 
to be identically 1/2; thus, the partitions of the graph into cycles correspond to 
the ways of choosing r on the interior faces such that whenever fi and f 2 are 
neighbors, |t(/i) — r(/2)| = 1/2; and if / shares an edge with the unbounded 
face, / differs by plus or minus 1/2 from the boundary value of that edge. In 
our grid graph illustrations in the next section, we multiply r by two so that we 
can deal with integer boundary values and increments of one. 

4 Computing the Hamiltonicity of Grid Graphs 

4.1 Fixed Faces and Grid Graphs 

In this section, we let G be a grid graph, i.e., an induced subgraph of Z 2 , and 
take <f> to be identically 2. Thus, S$ consists of the spanning subgraphs H of G 
that partition G into some number of cycles. We assume G is a polygon graph 
and define a hole in the graph G to be a bounded face of G with more than four 
edges. (All two-connected grid graphs are polygon graphs, and graphs that fail 
to be two-connected cannot be Hamiltonian, so we lose no generality with our 
assumption.) 

The first question we ask is, what are the connected components of if 
we only allow ^-transformations on the faces of G that are squares? Suppose 
H\ and H 2 are two partitions of G into cycles. Then uijjx — uj h 2 is closed and 
hence exact. Choose r so that dr = loh 1 — w/f 3 and r(/o) = 0, where /o is the 
unbounded face. A clear necessary condition for Hi and H2 to be connectable 
by ^-transformations on squares is that r be equal to zero on all of the holes of 
G. We claim that this is also sufficient. 

Theorem 8 Using the above notation, Hi and H 2 are in the same component 
of S(j) if and only if r(f) = for all holes f. 

Suppose without loss of generality that r(/) > for some square /, and let 
F + be the set of all faces of G on which r achieves its maximal value; note that 
by assumption, all of these faces are squares, not holes. 

Recall that if two faces fi and f 2 share an edge e, (/1 on the left with e 
is oriented black to white), then cjjj(e) can take on two values depending on 
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whether e G H. We say fi is higher (lower) than f\ (with respect to H), or 
/i -< /2, if wjf(e) takes on the larger (smaller) of these two values. We will also 
say that / is a local maximum (minimum) if it is higher (lower) than all of its 
neighbors. Clearly, one can apply a Z transformation to H at / if and only if / 
is a local maximum or minimum with respect to H . 

We will now show that we can always find a face f G F + such that / is 
higher than all of its neighbors, and hence, by applying a Z-transformation to 
Hi at /o, we can decrease |t| by one. First, define h(f) to be the number of faces 
adjacent to / that are higher than / minus the number that are lower. Because 
the contributions from edges between squares of F + cancel out, X)/gf + M/) 
is the number of edges of squares in F + that border on faces not in F + . In 
particular, the sum is positive, and the average value of h(f) is greater than 
zero. And this implies that there exists at least some fo in F+ which is higher 
than at least three of its neighbors — say, those on the south, west, and east. 
If fo is also higher than its neighbor on the north, we are done. If not, then 
let fi, be the face immediately north of fo', because it is higher than fo, it is 
necessarily in F+ and hence necessarily a square. 

Like fo, fi is higher than its neighbor to the south. We claim that f\ is also 
higher than its neighbors on the east and west. If, say, the western neighbor were 
higher, there would be an increasing sequence f w ,fo,f,f of faces all incident 
to a single vertex v; this would require either that all three of edges dividing 
these faces be in H or all three not in H. This is impossible since <j>(v) = 2. 

Continuing, we choose /b, /3, • • ■ until we find a face fk in F + that is higher 
than all of its neighbors. Because the sequence /i,/2, • • • must leave F+ even- 
tually, we know we find such an fk- □ 

Using the terminology of the previous section, this theorem is equivalent to 
the statement that are no face clusters of G that consist entirely of squares. 

Thus, if G has m holes, then when we only allow Z-transformations on 
squares, the connected components of S<j, are indexed by the integer lattice 
points in a convex polyhedron of dimension at most m. How many of these 
components are there? Well, if we take F to be the set of square faces of G, 
then there is a path from any hole to the boundary face or some other hole 
that passes through at most y/\F\ faces. We say two holes (or a hole and the 
boundary face) are "close" if they are adjacent or if there is a path from one 
to the other that passes through at most y/\F\ squares. One can easily show 
that the set of holes and the unbounded face are connected under the closeness 
relation. The number of values that th can take on any given hole once the 
value at a close hole is determined is bounded by 2-y/|F|, and it follows that 
there are at most 2 m \F\ m / 2 connected components of S^. 

In the remaining subsections, we will only deal with ^-transformations on 
squares, not holes, and we use only these transformations to describe adjacencies 
between members of S^. Our approach to determining the Hamiltonicity of G 
will be to search through each connected component of for a Hamiltonian 
path. Let p(H) be the number of cycles into which H partitions G. In this 
section, we will write H ~ Ho when H and Ho are in the same component of 
S,/,. (Equivalently, H ~ H when ujh and lJh are homologically equivalent in 
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the plane minus the holes.) 

In coming sections, we assume that we begin with an Hq G and will show 
that in 0(|F| 2 ) steps, we can find the H ~ H for which p(H) is minimal. It will 
follow that in 0(\F\ m / 2+2 ) steps, we can find an H for which p(H) is minimal 
overall — in particular, we will know whether G is Hamiltonian. 

4.2 Height Functions for Illustrations 

If G doesn't have any holes, then instead of considering the closed 1-chain ujh — ol 
and the r that has this as its boundary, it will be more instructive for us to 
consider the 1-chain Wh = 2w# — ft, where (3 is the sum of all edges in G oriented 
from black to white. Thus, if e is oriented black to white, u(e) = 1 if e is in H 
and —1 otherwise. 

Since du>(v) is not necessarily zero when v is on the unbounded face, we must 
adopt the "boundary value" approach described earlier; that is, we compute the 
boundary value for all the edges of the unbounded face. (If two boundary edges 
share a vertex v, the difference between the value of any th on the two edges is 
determined from the fact that <j>{v) = 2.) Now the set of allowable th is the set 
of functions on the other faces of G whose value on each face differs from the 
value of its neighboring faces (or boundary edges) by 1 or — 1. The mnemonic 
to keep in mind is that as we move from square to square with a black vertex 
on our left, tjj goes up if e is not in H and down if e is in H . 

When G has holes, the problem of producing tjj with ZJh as its boundary 
is analogous to the problem of integrating a complex function in a multiply 
connected region. Even in this case, however, we can define th on polyomino 
subgraphs of G and treat it as a locally definable multivalued potential function. 
We will label squares of most of the examples to come with the integer values 
of a locally defined th ■ This illustrative approach will make it easy to see when 
squares are locally maximal or minimal and when they are higher or lower than 
their neighbors. Note also that applying a Z-transformation at a square / 
corresponds to changing r#(/) by 2 or —2. 

4.3 Joining Cycles with Z-transformatons 

Let B{H) be the set of boundary faces, that is, the faces of G which contain 
vertices from two or more distinct cycles. (These are not to be confused with the 
primarily illustrative notion of boundary edges.) Now, if / is a locally maximal 
or minimal square, let H' be the graph obtained from H by applying a Z- 
transformation at /. Then if / G B(H), applying the Z-transformation at / 
joins the cycles on opposite sides of / together; in this case p(H') = p{H) — 1 
and / ^ B(H'). Similarly, if / B(H'), applying at a Z-transformation at / 
divides the cycle with vertices on / into two cycles; in this case, / £ B(H') and 
p(H') = p{H) + 1. The reader may check this fact on the local maxima and 
minima in the following example. 

Since we seek the H for which p(H) is minimal, the obvious first step is to 
hill climb until we either find an H for which p(H) — 1 (and we are finished), or 
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Figure 1: Applying a ^-transformation at a boundary (non-boundary) square 
decreases (increases) the number of disjoint cycles by one. 



we find an H such that B(H) contains no locally maximal or minimal squares. 
What can we say about H in the latter case? 

Suppose B(H) contains no locally maximal or minimal squares, and consider 
a face / £ B{H). We say a vertex v on / is isolated if neither of the vertices 
of / adjacent to v is in the same cycle as v. If two isolated vertices v\ and v 2 
that are incident to / are adjacent, then they belong to different cycles, and 
the face which shares the edge (vi,v 2 ) with / is necessarily a locally maximal 
or minimal square in B(H), a contradiction. If / contains vertices from three 
or four cycles, it must contain an adjacent pair of isolated vertices; hence, / is 
on the boundary between exactly two cycles. 

Suppose / contains one pair of vertices from each cycle. If v\ and v 2 are 
from the same cycle, then V\ and v 2 must be adjacent vertices. (Otherwise, they 
would be isolated.) Similarly, the other two vertices V3 and V4 are adjacent. If 
both of the edges (vi,v 2 ) and (^3,^4) were contained in H, f would be locally 
maximal or minimal. If neither were contained in /, the two faces adjacent to 
/ containing vertices from both cycles would be locally maximal or minimal. 
Hence, H contains exactly one of (v\,v 2 ) and (^3,^4). Squares of this type are 
called critical boundary squares. Inspection shows that they are adjacent to 
exactly two boundary squares (or holes), and that they are either higher than 
both or lower than both boundary neighbors. Hence, we say they are maximal 
or minimal along the boundary. 

Similarly, suppose / has one vertex v\ from one cycle and three {v 2 , V3, 
and V4) from the other, with v 2 and v 4 adjacent to v\. Since boundary squares 
adjacent to / cannot contain edges from both cycles without being maximal 
or minimal, one easily checks that (^2,^3) and (^3,^4) must be in H. Squares 
of this type are called corner boundary squares, and they are not maximal or 
minimal along the boundary. We summarize this information in the following 
lemma. 
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Lemma 1 If B{H) contains no maximal or minimal squares, then every square 
in B(H) is either a critical boundary square or a corner boundary square, as 
illustrated in the following diagram. 



2 3 2 4 3 2 

Figure 2: Critical and corner boundary squares. 



The examples show two graphs with multiple cycles for which \S$\ = 1 , and 
hence joining the cycles is impossible. This is because each square in these 
graphs is on an increasing path from one boundary edge to another; intuitively, 
the height functions are pulled taut by the boundary values, and there is not 
enough slack to join the separate cycles together. 




Figure 3: |S | = 1, \B(H)\ = 6, and p(H) = 2 



4.4 Height Values on Boundaries 

We begin this section with a simple observation about cycle boundaries: 
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Figure 4: |S* | = 1, \B(H)\ = 32, and p(H) = 3 



Lemma 2 If the function th> has the same value as t h on every boundary 
face in B(H), two vertices that are on separate cycles in H are necessarily on 
separate cycles in H' . 

Otherwise, there would have to be an edge in H' connecting a pair of vertices 
which belonged to different cycles in H. This cannot be, because the faces f\ 
and fi on opposite sides of that edge are in B(H), but the difference in values 
on the two faces is not the same for m and th> ■ D 
Another simple observation is the following: 

Lemma 3 If th> has the same value as th on every hole and every critical 
boundary square, it necessarily has the same value on every face in Bh- 

To see this, note that every corner square is adjacent to two boundary squares 
or holes with which it doesn't share an edge in H ; one of these is higher and one 
is lower; if these are also corner squares, they are in turn adjacent to respec- 
tively higher or lower boundary squares or holes. The sequence of increasingly 
higher corner squares, will step alternatively in two directions (say north and 
east) while the downward sequence alternates between the other two (say south 
and west). Clearly neither sequence ever returns to a previous square in the 
path; hence the upwards (downward) sequence must eventually reach a locally 
maximal (minimal) boundary square or a hole. Thus, each corner boundary 
square is on an increasing or decreasing path between holes or squares that are 
minimal or maximal along the boundary. □ 

From these two lemmas, we deduce that if we hope to apply a sequence of 
^-transformations to H to produce an H' such that p(H') < p{H), we will at 
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sonic point have to apply a ^-transformation that lowers (raises) at least one of 
the critical boundary values of H which is currently maximal (minimal) along 
the boundary. The next section describes a way to do just that. 

4.5 Maximal and Minimal Rows 

Suppose that fo is a critical boundary square and that its edge e containing 
vertices in one cycle is in H and its edge e' containing vertices in the other 
cycle is not. Suppose that e has / on its left when oriented from black to white. 
(The other case — as should be understood thought this section — is symmetric.) 
Then / is higher than all of its neighbor faces except for f\, the other face 
incident to e'. Let fi, f 2 , ■ ■ ■ be the sequence of faces extending away from / in 
the direction of the edge e' . Let r be the smallest integer such either f r is no 
longer a square or f r is lower than f r -i- Clearly, r > 2. 

If f r is a boundary square or not a square at all and f r is higher than / r _i, 
then we say that / is a stuck critical boundary square. In this case, there is no 
way we can lower fo without first lowering f r , which is either another boundary 
square or a hole. 

Otherwise, the sequence fi, ■ ■ ■ , f r -i consists of squares all of whose vertices 
belong to the same cycle. Also, /o,/i,... ,/r-i is a maximal square row of 
length r, that is, a sequence of r consecutive squares in a row or column such 
that fi -< for < i < r — 2 and each square in the sequence is higher than all 
of its neighbors not in the sequence. We can lower the maximal square row by 
decreasing the value of th by two on all elements in the row — or equivalently 
by applying our ^-transformations successively to / r -i, /r-2, • ■ • ,fi,fo — to 
produce a new subgraph H'. 

The effect of this lowering is different depending on whether r is even or odd, 
as shown in the following figures. If r is odd, this has the effect of joining the 
cycles separated by fo; in this case p(H') = p(H) — 1, and B(H') is a proper 
subset of B{H). Moreover, th and th> agree an all elements of B(H'). 

If r is even, the process also joins the cycles on opposite sides of /. However, 
by assumption, the maximal square at the top of the maximal square row is 
not on a boundary. (If it were, we could easily join two cycles by lowering 
simply that square.) Hence, lowering that square creates an additional cycle. 
In this case, p(H') = P(H), and neither B(H') nor B{H) is a proper subset 
of the other. However, it is still the case that t(H') and t(H) agree on the 
intersection of B(H) and B(H'). 

Does lowering an even length maximal row get us anywhere? Maybe. It 
doesn't change p(H), but it may be the case that it creates a new boundary 
square that is the beginning of an odd length maximal or minimal row (possibly 
of length one). 

4.6 Bridges 

In the scenario of the last section, suppose r is even and /o, • • • , /r-i is a max- 
imal vertex row beginning at a boundary square fo- (By construction, / r _i is 



18 



7 7 




] 1 



Figure 5: When a boundary square is the first square in an odd length maximal 
vertex row, lowering that row joins the cycles on opposite sides of the square. 
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i L 



Figure 6: When a boundary square is the first square in an even length maximal 
vertex row, and the last square in the row is not on a boundary, lowering the 
row joins the cycles on opposite sides of the first square and divides the cycle 
with vertices on the last square into two pieces. 
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not on a boundary.) Then when we lower the vertex t> r -i, this partitions the 
cycle with vertices on / r _i into two cycles. At this point, / r _i is a minimal 
square on a boundary between two cycles and /o, • • ■ , f r -2 is an odd length 
maximal square row starting on a boundary between another pair of cycles. I 
refer to a row of squares /o, ■ • ■ , f r -i of this type as a bridge; that is, I say there 
is a bridge of length r from / to / r _i if / and / r _i are boundary squares 
and the squares /o, fi, ■ ■ ■ , are lined up in a row so that /o, • • • , /r- 2 is an 
odd length maximal row of squares and f r ~i is a minimal vertex row (of length 
one). If I raise f r -i, thus joining the cycles on opposite sides of that boundary 
point, then /o, /1, • ■ • , f r -i is no longer a maximal vertex row. If I lower the 
maximal vertex row /o, /1, . . . , /r-2, thus joining cycles on opposite sides of the 
boundary point /o, then is no longer minimal. Thus, I can think of a bridge 
as an option to join cycles at one square or another but not both. 
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Figure 7: Bridges of length six, four, and two. 



After constructing a bridge, I refer to the squares / and f r _i as bridge 
end squares and the squares /1, • • • , f r -i as inner bridge squares. Creating and 
using these bridges will be central to our algorithm. 

4.7 The Algorithm 

This section describes a method for reducing the number of cycles of H by one; 
we will show that this method fails only when H has the minimum number of 
cycles among members of its connected component of S<j, . 

We assume H contains at least two disjoint cycles, and by the boundary 
between two cycles, we will mean the set of all faces that are incident to both of 
those two cycles. We begin the algorithm with B\, . . . ,Bk, all of the nontrivial 
boundaries between pairs of cycles. While wc arc ultimately searching for an 
odd length maximal or minimal row that will allow us to get rid of one of these 
boundaries, during the course of this algorithm, we may temporarily create new 
boundaries Bk+\, ■ ■ ■ , B m . When a given Bi is created, it represents the set of 
squares between exactly one pair of cycles; as the algorithm progresses, these 
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may later be separated into subcycles, but that will not concern us. We need 
only remember that each Bi was a boundary between exactly two cycles when it 
was created. Also, each time I create a new boundary, I will also create a bridge, 
as described in the previous section, from a square on a previously existing cycle 
boundary to a square on the new boundary. 

The algorithm works as follows. We start successively searching all of the 
critical boundary squares /o of H. If /o is stuck, or / is a square on a previously 
created bridge, we do nothing. Otherwise, /o must be the beginning of a maxi- 
mal (minimal) row of squares of length r. If r is odd, then we lower / r _i, thus 
producing a new boundary containing / r _i and a new bridge from /o to / r _i. 
We say that the new cycle boundary (the one containing / r _i) is subordinate 
to the cycle boundary containing / . We continue this process until one of two 
things happens. 

1. We find a critical boundary square /o (perhaps on one of the new cycle 
boundaries) that is the beginning of a maximal or minimal row of even 
length. 

2. Every critical cycle boundary square is either a stuck critical boundary 
square or a bridge square. 

Ultimately, we'd like to show that if the first case occurs, we can lower or 
raise vertex squares in a way that decreases p(H) to one below its original value. 
To make this work (and simplify our proofs), we will need two additional rules 
about the order in which we select our critical squares. Since we'd prefer not to 
deal with squares that are boundaries between more than two cycles, if we ever 
produce a new boundary that contains an / on the boundary between three 
cycles, then the next square we check in our algorithm is /', the square adjacent 
to / across an edge connecting two isolated vertices of /. We saw earlier that 
this square is necessarily maximal or minimal; one can easily check that it does 
not belong to any existing bridges. This fact means that we will not have to 
deal with boundary squares between three or more cycles at any other point in 
our algorithm. Except for this one case, each time a new Bi is created, if one of 
its squares is adjacent to a square in another B i: the two are separated by an 
edge in H. There is one other rule that we will add when it is motivated later 
on. 

During this process, we will have organized our new boundaries into trees 
- under the subordination relation - rooted at the original B\ , . . . , Bk ■ Each 
time we add a new boundary, we add it as a leaf to an existing node. We will 
now show that once we find a cycle boundary square on an even length maximal 
(minimal) row, we can join cycles in such a way as to eliminate all of the new 
boundaries -Bj+i, • • • , B m we have created and the boundary B t which is the 
root of the boundary tree containing B k , thus producing a partition of G into 
cycles with one fewer cycles than our original partition. 

We use the following diagram to give an overview of this process. The roots 
of the trees B\, . . . ,B k represent the original boundaries between pairs of cycles. 
The -Bfc+i, ■ ■ • ,B m represent new boundaries that are erected (in that order) 
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during the algorithm. Recall that a bridge is made up of a minimal and a 
maximal odd length row; and we can lower or raise exactly one of these rows to 
join cycles across the boundary at one of the bridge's endpoints. Each edge in 
the diagram represents a bridge across one of the two boundaries it is connected 
to — except for the loose edge off of B m which represents a maximal (minimal) 
row of squares of odd length starting at / that we will call a half bridge. 



Figure 8: Using bridges and one half bridge to eliminate boundaries. 

The arrows on the bridges point to the boundary across which we will use 
the bridge to join cycles. As shown, at each bridge we will join cycles across the 
subordinate boundary except for the bridges on the path from a root boundary 
to the half bridge; as shown in the diagram, we will then have joined cycles 
across all of the newly created boundaries and the root boundary Bj. In order 
to make this precise, we will need the following lemma, which we will prove 
later. Let Mj be the even length maximal (minimal) row of squares which is 
lowered (raised) in order to join cycles across the boundary Bi. 

Lemma 4 As sets of squares, all of the Mi used by our algorithm are disjoint. 
Furthermore, no square in any one of the Mi is adjacent to any square in another 



The lemma shows that lowering or raising one of the Mj will not affect the 
minimality or maximality of any of the other Mi. For the remainder of the 
proof, it is useful to imagine that we join cycles across the boundaries one at a 
time in the reverse order of the order in which we created them, that is, starting 
with B m , progressing to Bk+i- 

When B m was created, it separated a particular cycle into two pieces; by 
joining the two cycles on opposite sides of B m , we reduce B(H) to exactly 
what it was before B m was created. We do the same thing by joining cycles 
on opposite sides of B m _i, B m -2, • • • B/.+1- To see this, observe that the bridge 




Mi. 
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that we use to join cycles across a given Bi is either the bridge used to create 
that Bi or a bridge produced later in the algorithm, it follows that the bridge 
was created while all B n , n < i were in place — so by our construction, using 
it to join cycles across Bi cannot alter the value of th on any other of the B n 
and thus will not alter B(H) other than by eliminating the boundary squares 
of Bi. 

Hence, after joining cycles across Bk+i, B(H) is the same as it originally 
was. We then join cycles across Bi to reduce the number of cycles in £ by one. 
□ 

Proof of Lemma The proof of this lemma is slightly tedious, but not 
difficult. It involves checking a number of cases. First, we will show that if a 
bridge /o, • ■ • , f r -i is oriented in the vertical direction, none of the other bridges 
or half bridge can contain any of the squares adjacent to /i, . . . , f r —i on the 
left or right. (Of course, the horizontal case is equivalent.) 

We prove this by showing that each time we construct a new bridge or half 
bridge, it cannot contain squares from any other bridge, and moreover, that if 
an element of one bridge is adjacent to that of another bridge or half bridge, 
the elements belong to maximal (or minimal) rows that will not both be used 
(i.e., lowered or raised) by our algorithm. 

So now we look at what happens when we construct a new bridge or half 
bridge. The algorithm requires us to find a cycle boundary square f that is 
part of a maximal row /o, /i, ■ ■ • , ,f r -i, and we then either lower f r -i (in the 
even case) to produce a new boundary and a bridge or the entire row (in the 
odd case). Now, recall that before we do this, none of the /i, . . . , f r -i is on a 
boundary, and it is thus possible to lower f r -\ through /i without changing the 
value of th on any boundary square. 

Now, suppose g, g\, . . . , gu-i is an existing bridge oriented from top to bot- 
tom as in the following diagram. As can be seen from the diagrams, the edges 
separating any squares immediately to the right or left of the g,g±, . . . , gu-i are 
alternatively contained and not contained in H, according to the parity of their 
endpoints along the bridge. It follows that these squares are all on directed 
paths from one boundary square or hole to another, and we thus cannot apply 
a Z-transformation to any of these vertices without first changing the value of 
th on a boundary square. Thus, none of these will be among the fi, ■ ■ ■ , f r -i- 

Hence, the only way our new bridge or half bridge can contain a square 
adjacent to one of the g, . . . , gu-i on the left or the right is if that square is the 
starting square, fa. Then /o would have to be adjacent to cither g or gk-i, since 
otherwise it would be in the middle of a directed path from one cycle boundary 
square to another, and hence, if critical, it would have to be stuck. The new 
bridge then has to be oriented vertically in this case since it begins on a boundary 
square with another boundary square as a horizontal neighbor. The two bridges 
are adjacent only at that one square. It is easy to check from the diagram that 
in this case / and g must be on the same Bk , because they necessarily separate 
the same pair of cycles. (This follows from the fact that there are no triple 
boundary squares; that case is dealt with separately.) However, only one of 
these squares (say /) can be part of Mk] the other bridge must be used to join 
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cycles across a different boundary, say Bi, and g will not be part of Mj. 

Now, if the new bridge or half bridge is oriented horizontally, we have shown 
that none of fo, ■ ■ ■ fr-i can be adjacent to a vertical bridge on the left or right, 
and hence none of /o, . . . , f r —i can overlap with a square in the vertical bridge. 
Similarly, if the new bridge is oriented vertically, fo cannot (by construction) be 
one of g, . . . , gu-i, and neither of g and gu-i can be among the /o, • ■ • , f m -\ 
(by construction); hence, the squares in the bridges cannot overlap in this case 
either. 

We still have one case to deal with, namely, that the new bridge /i, . . . , f m -\ 
may contain either the vertex gk immediately above gu-i or the vertex g_i 
immediately below g. First we rule out the horizontal case. One can easily check 
that both before and after a horizontal bridge is created, all of the squares that 
belong to that bridge contain at least one pair of vertices from which vertical 
edges extend away from the square in either direction. If either g_i or gk were to 
be part of a bridge oriented left to right (or any part of a half bridge except the 
last square) , there would have to first be edges extending upward and downward 
out of those squares, which cannot happen since those edges would have to cross 
boundaries; these squares cannot be the last element of a half bridge oriented 
left to right because they border an incident edge in H either above or below. 

Since the new bridge or half bridge cannot overlap with the g, . . . , gk-i, the 
only remaining possibilities are that one of gk or g_\ is a starting point or ending 
point for a vertical bridge or half bridge. Assume that vertex is g-\. (The other 
case is identical.) 

Now, the bridge creation cannot raise or lower g-i, since this would neces- 
sarily change it to a vertex with edges on the right and left and, if a bridge end, 
it would thus have to be the end of bridge oriented horizontally. It follows that 
<7_i has to be either the starting square of a bridge, the starting square of a half 
bridge or the ending square of a half bridge. 

For these final three cases, we will change our strategy. Instead of showing 
that these three situations can't possibly arise (in fact, they can), we will show 
that whenever our algorithm does attempt to include g_i in one of these roles, 
there is a maximal or minimal square to the left or right of g that we can use as 
a half bridge instead of creating the bridge or half bridge containing (This 
is the second rule alluded to earlier.) 

First of all, in all three of these cases, the incident edges to the left and right 
of g-i cannot be in H. It follows that the edges above the squares left and right 
of g-\ must be in !K, as seen in the following figure. 

Since g_i was to be part of a new bridge or half bridge, it follows that neither 
a_i nor 6_i may be an inner bridge square. Thus, if there is a bridge containing 
a (resp. 6), it must be a bridge of length one from a to a_i (resp. b to 6_i). 
Also, we cannot have bridges at both of these places, since they would then have 
to join the same pair of boundaries, and this never happens in our algorithm. 
(If g_i were a triple boundary square, then a_i and &_i would not necessarily 
separate the same pair of cycles; however, by assumption, this is not the case.) 

Now, suppose there is a bridge from a to a_i and not from b to b-\. Then 
we take the square b to be our half bridge. Let Bk be the boundary containing 
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Figure 9: A bridge connects g and gk-i when the algorithm attempts to include 
g-i in a bridge or half bridge. The squares a_i and &_i may actually be parts 
of holes, but the edges between a and a_i and between b and b-\ are necessarily 
present. 



b — then consists of the single square b. Now, we must show that b cannot 
be adjacent to squares in any of the other M,. We have already dealt with the 
case of the squares on the left and right of b. 

Now, suppose that b\ is a bridge square on a boundary £?/. This implies 
that the bridge starting at g must be of length two, since otherwise, b\ would 
be adjacent to an inner bridge point. Since we are joining cycles across Bk with 
6, we must use the bridge containing g and g\ to join cycles (at wi) across £>;. 
Thus, the bridge starting at &i will not be used to join cycles across Bi, and 
hence 61 will not be in any Mi. 

Similarly, suppose that 6_i is a bridge point on a boundary Since we 
are joining cycles across Bk with 6, we must use the bridge containing a and 
a_i to join cycles (at a_i) across B\. Thus, the bridge starting at 6„i will not 
be used to join cycles across Bi, and hence 6_i will not be in any Mi. These 
figures illustrates some of these possibilities. 




Figure 10: Examples of other bridges that may be present. 
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Now, the remaining case is when there is neither a bridge from o to a_i nor 
from 6 to 6_i. In this case, if there is no bridge starting at a_i (rcsp. 6-1), we 
can take our half bridge vertex to be a (resp. b), and then the proof that this 
square is not adjacent to any other Mj is the same as in the previous case. 

However, suppose there arc bridges beginning at both a_i and 6_i, and these 
two squares are contained on the boundary B\. Then at most one of the two 
can be part of M;, and we can determine which one that is from the boundary 
tree structure using the fact that our half bridge will be at Bk- Suppose a_i is 
a vertex that is not part of Mi (and hence not part of any Mi). Then let a be 
our final half bridge vertex, and the proof that a is not adjacent to any vertex 
in one of the Mj is then the same as in the previous cases. 

As a final note, we have already shown that the square / we use as our final 
half bridge after producing a triple boundary square is not a bridge square. 
Inspection shows that the squares adjacent to / across edges contained in H 
cannot be bridge squares since they cannot be critical boundary squares. □ 

4.8 When Algorithm Gets Stuck, Number of Cycles is 
Minimal 

Suppose that in the above algorithm, we started with m cycles, then added r 
new cycles and r new bridges, and then got stuck. That is, we produced an H 
with the following properties: 

1. H divides G into m + r cycles, and there are exactly r non-overlapping 
bridges designated. 

2. A square of one bridge may be adjacent to a square of another only if 
the two bridges are oriented vertically (horizontally) and the two squares 
are both bridge ends; no two bridges are ever adjacent at more than one 
square. 

3. Every critical boundary square is either a stuck square or a square on a 
bridge. No boundary squares have vertices from more than two cycles. 

Now, for each bridge /, /1, . . . , / m -i, we will refer to the squares f m -2, fm-i 
as a fixed pair. The fixed pair looks like the length-two bridge shown in a 
previous diagram; that is, the two squares are separated by an edge in H, 
the higher of the two is a local maximum, and the lower is a local minimum. 
Intuitively, these pairs occur at places where there is slack in the height function 
that could be used to join cycles together; however, we cannot make use of that 
slack because we are holding the edge between those squares to be fixed. 

Now, our system satisfies: 

1. H divides G into m + r cycles, and there are exactly r designated fixed 
pairs which do not overlap. 

2. Whenever a square of one fixed pair is adjacent to a square of another, the 
two fixed pairs are both oriented vertically (horizontally). No two fixed 
pairs are ever adjacent at more than one square. 



26 



3. Every critical cycle boundary square is either a stuck square or a member 
of a maximal (minimal) row of squares, whose latter endpoint is a member 
of a fixed pair. 

We would like to prove that any other H' in the same connected component 
of S,p has at least m distinct cycles. We will do this by showing that we can 
in a number of steps alter H in such a way that each step will preserve all the 
properties listed above, and so that when we are done we will have produced a 
new H' with potential th> , such that th> agrees with th on all the squares on 
all but r of the cycle boundaries of if'. Joining the cycles across a particular 
one of these boundaries can reduce the number of cycles in if by at most one, 
so we will conclude that if' has at least (m + r) — r = m separate cycles. 

Suppose that t' h and th are not in complete agreement, and let 

F+ = {f\Mf)>TH>{f)} 
F-={f\T H (f)<T H ,(f)} 

Now, suppose that / G F + (the F_ case being symmetrical) and / is maximal 
along the boundary. That is, / is a critical square and / is higher than its two 
boundary square neighbors. Then either / is maximal overall (in which case 
it must be a bridge square), or there is some /i higher than /. In the latter 
case, if /, ,f m -i is a maximal vertex row starting at /, then each of 
/i, . . . , f m -i must be in V+. Since th and t' h agree on holes and the unbounded 
face, /i, . . . , f m -i is a maximal row of squares, and hence / m _i is a member 
of a fixed pair. Since neither square to the left or right of f m -\ is a local 
maximum or minimum, this fixed pair must be oriented the same direction as 
the /i, . . . ,/m-i- 

Now, as the diagram shows, we will remove the fixed pair containing / TO _i, 
lower the vertices /i, . . . , f m -i, and add /, fi as a new fixed pair. One easily 
checks that this new system satisfies the same three properties as the old one. 
The first property is obvious. The second is also fairly obvious from the dia- 
gram once one notes that a vertically (horizontally) oriented fixed pair is never 
incident to any vertical (horizontal) edges. This fact rules out a horizontally 
oriented fixed pair adjacent to the new fixed pair. A vertically oriented pair 
cannot be adjacent to the new fixed pair at two places on the left or right be- 
cause one of the squares on each of the left and right is necessarily not maximal 
or minimal. The third property also follows from inspection once we notice that 
each of the squares to the left or the right of the central column in the diagram 
is on a directed path between holes and/or boundary squares that were already 
critical (hence, newly created critical boundary squares on the sides don't pose 
a problem); also, critical squares along the row containing the new fixed pair 
that were previously part of maximal or minimal rows with an endpoint in the 
old fixed pair will now be part of maximal or minimal rows with an endpoint 
on the new fixed pair. 

We repeat this process for all other critical boundary points in F + and 
similarly in F- until every critical boundary square is either 
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Figure 11: When a critical boundary square is part of an odd length maximal 
row of length at least three and its other end square is the first of a "fixed pair," 
we can lower all the squares except the critical one and move the fixed square 
from the top to the bottom. 



1. A square on which th and th> are equal. 

2. A member of a fixed pair in F + that is maximal. 

3. A member of a fixed pair in F_ that is minimal. 

4. A square in F + that is minimal along the boundary. 

5. A square in F_ that is maximal along the boundary. 

If fi and are members of a fixed pair, and f\ is in F + and is maximal, then 
although fi may not be in F + , it cannot be in F_ because it is adjacent to f\. 
Similarly, if f\ is maximal in F + , fi cannot be minimal in F- . It follows that 
the number of boundary squares in F+ that are maximal along the boundary 
plus the number of critical points in F_ that are maximal along the boundary 
is at most r, the fixed square pairs. 

Next, let Bi be a connected component of the set of boundary squares sepa- 
rating some single pair of cycles. Because there are no squares separating three 
distinct cycles, one easily checks that every boundary square must be adjacent 
to exactly two faces (which may or may not be squares) on the same boundary, 
so that this set of squares can be written as a sequence / , . . . , fk such that 
each face is adjacent to its successor and either f and fk are adjacent to holes 
or they are adjacent to each other. In both cases, the edges between adjacent 
squares are not contained in H . 

Next, we claim that if Bi has a square in F + (F_ , it must contain a square in 
F + (F_ ) , that is maximal (minimal) along the boundary. If fi is the first square 
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in F + (necessarily fi >- f%-\), then the first square after fi which is higher than 
its successor is a local maximum; there must be such a place, because th and 
th> agree on the endpoints. 

If Bi contains only squares but contains one or more squares in F + and one 
or more not in F + , then we can apply the same argument to the sequences 
of squares contained in F + . If all of the squares in Bi arc in F + , then it is 
sufficient to note that every closed sequence of boundary squares must contain 
both maxima and minima. To see this, let / be any vertex on the boundary 
with a neighbor /i such that fi >~ f- Then either fi -< -< -< ...-</ or 
there is some vertex that is a maximum along the boundary. However, if in fact, 
fi ~> fi ~< /3 ■ • ■ then every square is a corner square, and the sequence must 
progress diagonally and never return to /. 

Thus, every connected component of the cycle boundary structure that does 
not contain a vertex in F+ that is maximal along the boundary or a vertex in 
F- that is minimal along the boundary must have th and th> equal everywhere. 
Since there are at most r vertices of this type, there are at most r connected 
components of the cycle boundary structure on which m and th< disagree. 
Joining cycles across each of those r components reduces the number of cycles 
by at most one, so m is indeed the minimum number of cycles for the connected 
component of S$ containing H. 

4.9 Runtime Analysis and Implementation 

The most naive implementation gives a runtime of at most (9(|F| 3 ) in the worst 
case. This assumes we start out with 0(|F|) separate cycles. Each time we try 
to join a pair of those cycles, we produce 0(|.F|) bridges. Each time we produce 
a bridge, we spend 0(\F\) time updating our list of critical boundary squares. 

To produce the bound 0(|F| 2 ), we still assume we start out with 0(|F|) 
separate cycles, but we claim that given H , we can always reduce p{H) by one 
(if possible) in 0(|F|) steps, because it is not really necessary to spend 0(|F|) 
steps updating our list of critical boundary squares each time we produce a 
bridge. 

We show this first in the polyomino case. First, we compute B(H) and 
put all the critical boundary squares in a list in 0(|F|) steps. Next, each time 
we lower or raise a square / to produce a bridge, we create a new boundary 
containing /. Now, if we knew that this new boundary contained a locally 
maximal or minimal square, we could take that square as our new half bridge; 
regardless of where it was, we could recompute B(H), find that square, and 
complete this phase of the algorithm in 0(|F|) steps. But for now, we will 
assume the boundary contains no locally maximal or minimal squares and try 
to reproduce it. 

That is, we attempt to form paths of non-minimal, non-maximal squares 
starting at / and extending in either direction such that the paths do not pass 
through any edges in H; we stop when either the paths meet each other (and 
we have a cyclic path) or they reach a previously existing boundary square 
or the unbounded face. The key observation is that under these assumptions, 
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there is only one way to create the path. To see this, suppose was one 
tentative boundary square, and fk+i was my next choice, immediately north 
of fk- Assuming fk+i is not a dead end, H contains either zero, one, or two 
edges of /fc+i; if it contains two edges, there is only one choice for fk+2- If it 
contains one or zero edges, there are more choices; however, assuming fk+i is 
not locally maximal or minimal, quick inspection shows that all but one of these 
is necessarily a dead end or a local maximum or minimum. 

If my sequence gets stuck in a dead end, I conclude that I must have taken a 
wrong turn — and hence my boundary contains a local maximum or minimum, 
and I'm done. If the sequence runs into a previously existing boundary square, 
then cither I took a wrong turn (and there's a local maximum or minimum on 
the new boundary) or I've created a triple boundary square. If the sequence is 
forced to choose a minimal or maximal square, then cither that square is correct 
or there was another minimal or maximal square off of a wrong turn. In each 
case, there exists a half bridge somewhere, and I can find it and complete this 
phase of the algorithm in 0(\F\) steps. If the sequence meets itself or if both 
ends reach the unbounded face, then I conclude that I really have found the 
new cycle boundary squares, and I proceed with the algorithm. 

For another enhancement to this phase of algorithm, I note that each time I 
create a new bridge square, it remains a bridge square until I'm finished; I label 
these squares as bridge squares. Similarly, if I find a critical boundary square is 
stuck, I label the entire increasing row of squares between that square and the 
other boundary square or unbounded face edge as stuck. If any of those become 
critical boundary squares at later stages in the algorithm, I won't have to check 
them again. With this scheme, I check each square as part of a new boundary 
creation stage at most once and I check each square as part of a bridge creation 
or stuck critical boundary square determination step at most once. I conclude 
that this entire phase of the algorithm is 0(|F|). 

If G has holes, then when I create a new boundary, even if I don't create any 
triple boundary squares, the new boundary may still contain several connected 
path components, each connecting one boundary to another. This suggests that 
my previous technique of computing the boundary squares won't work, because 
when my path reaches a hole, it's not clear where to go next. 

As it turns out, we can use the same idea as in the polyomino case with 
minor modifications. Suppose there are k holes counting the unbounded face. 
(For the purpose of algorithm bounds, treat fcasa fixed constant.) Then we say 
two holes are adjacent if there is a path of squares from one to the other that 
doesn't pass through any boundary squares or cross any edges of H and such 
that no square along the path is maximal or minimal. As part of our 0(|F|) 
preparation for this phase of the algorithm, we compute B(H) and determine 
whether each pair of holes is adjacent, storing the path that connects them if 
there is one. (There is necessarily only one such path — if there were two, at least 
one would contain a boundary face.) We also make a table indicating when two 
such paths overlap. Observe that under this definition of adjacent, the graph 
on the holes is a disjoint union of trees. 

Now, after we lower or raise an / to produce a bridge, we try to guess the 
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new boundary much the same way as before. We produce paths of squares 
extending in either direction; if they close in on each other or they both reach 
the same hole, we have found our boundary. If one reaches a dead end or 
another boundary face, we conclude that there was a maximum or minimum 
somewhere and we find it. If the paths reach different holes, then we check to 
see if those two paths are on the same adjacency tree. If they are, we take the 
unique set of boundaries connecting them to be our new boundary. If they are 
not, we conclude that there is a maximum or minimum along our new boundary 
somewhere and we find it. 

5 Generalizing to Gridlike Graphs 

A careful examination of the steps in the above proofs reveal that the require- 
ment that G be a grid graph is much stronger than what we needed. In this 
section, we try to capture the "grid-like" properties that are actually necessary. 

Let G be a graph embedded in on a plane or n-holed torus such that all of 
its faces are simple polygons. We say a face / of G is a good square if it satisfies 
the following properties: 

1. Every vertex of / has degree at most 4. 

2. If two adjacent vertices V\ and v 2 of / both have degree four, then the 
face /' which shares the edge (^1,^2) with / is a (not necessarily good) 
square. 

3. The square / does not share more than one edge with any single square. 

A k-holed, n-holed-torus-embedded grid-like graph is a graph G embedded in 
an n-holed torus such that all of its faces are polygons and all but k of its faces 
are good squares. 




Figure 12: Planar grid-like graph with no holes. 

In this section, we say that H ~ H if u>h — ^h 1S exact and let r be the 
2-chain with this as its boundary (such that r(/o) = for a reference hole, 



31 



/o). We claim that the algorithm developed in the previous section extends to 
gridlikc graphs. That is, if we fix n and k to be constants, then given a fixed H 
in a fc-holed, n-holed-torus-embedded grid-like graph, we can find an H ~ H a 
for which p(H) is minimal in 0(|^| 3 ) steps. 

Both the algorithm and the proof of its effectiveness are the same as in 
the grid graph case with a few modifications. We could have introduced this 
notation at the outset, but that would have made the first read-through of the 
algorithm more confusing and only slightly more enlightening. So instead, we 
outline here the modifications required and present the details as exercises. 

1. Prove that if none of the good squares in B(H) is locally maximal or 
minimal, then all good squares arc cither corner boundary squares or 
critical boundary squares. (The first step is to use the definition of good 
square to show that if two isolated vertices v\ and v 2 on a good square / 
are adjacent, the face that shares the edge (v\,V2) with / is necessarily a 
square and hence maximal or minimal.) 

2. Construct maximal or minimal rows as follows. Let /o be a square that 
is higher than its neighbors on the south, east, and west, and lower than 
its neighbor to the north. Show that if it is possible to decrease Tn(fo) 
without changing th on any holes, then Tn(fo) must be part of a unique 
increasing path of good squares such that every square in the path is 
higher than every face outside the path. For i > 0, define the south edge 
of fi to be the single edge that borders fi-i, and define north, east, and 
west edges accordingly. Show that each fi in the sequence is higher than 
its neighbors to the east and west. 

3. In what ways can squares in a maximal or minimal row be adjacent to 
other squares in the same row? Answer this question with the following 
steps. Let i be the first integer such that fi shares an edge e with some 
earlier fj in the sequence (not counting its neighbor to the south). Note 
that e cannot be the east or west edge of fi because if it were, fi would 
have to be higher than fj, and fj is higher than all squares except for its 
neighbor to the north. Suppose e were the north edge of fi and the west 
(or, symmetrically, east) edge of fj. If /' is the neighbor of fj to the west, 
then /', /j, /i, is an increasing sequence of squares sharing the same 
vertex, and this also cannot happen. Hence, it must be the case that the 
north edge of fi is the south edge of fj, and thus j = 0, since already we 
know the southern neighboring face of the other fj. 

In this case, we refer to the sequence /o, ■ • • , fj as a cyclic maximal good 
square row. Since G is bipartite, it has even length (i.e., j is odd). It 
is impossible to change the value of th on these good squares with Z- 
transformations alone, but we can decrease the value of r on all the squares 
in this row at once. Show that this has no effect on the value of p(H). 

4. Show that if H ~ H , it is possible to move from H to H by successively 
applying ^-transformations and raising (lowering) cyclic minimal (maxi- 
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mal) good square rows. As before, this is done by finding an clement in 
F + (or F_) that is higher (lower) three squares. We saw earlier that this 
is possible provided F + borders at least one element not in F + . Show that 
this fails only when F + contains every face in F and there are no holes and 
every face is higher than exactly two squares. (Certain cycle partitions 
on C m x C n have this property. Show that in this exceptional case, S<p 
contains exactly two elements.) 

5. Show that all of the steps in the basic algorithm arc well-defined in this new 
setting. Adapt the proof of the lemma that says the maximal and minimal 
rows used to join vertices are disjoint and no element of one is adjacent to 
any element of another. One step requires showing that if /o, . . . , f m -i is 
a north-south oriented bridge, the faces adjacent to /i, . . . , f m -i on the 
left and right are on directed paths from one hole or member of Bh to 
another; in our new setting, we have to deal with the case that a single 
square / is adjacent to two or more of the /i, . . . , f m -i on the east or 
west. Still, we can let /g, /{, /V, . . . ,f m -i,f' m be the faces adjacent to 
/o, . . . , f m and show that even if the list contains repeats, it still describes 
an increasing path from one hole or boundary face to another. 

6. Using fixed pairs, show that when the algorithm fails, p{H) is minimal. 

6 Markov Chains for Hamiltonian Cycles 

Write = {H £ S r j > \p{H) = k}, i.e., the set of partitions of G into exactly 
k cycles. In this section, instead of simply using Z-transformations (raising or 
lowering single squares), we say two subgraphs of are adjacent if one can be 
obtained from the other by raising (lowering) a minimal (maximal) row. Let "K 
be a chosen homology class of S<t>. 

Theorem 9 Under the above definition of adjacency, (S^US^nJC is connected 
whenever both S^C\ H and «S| fl H are nonempty. 

To see this, we show that given an H' G S^, and an H e U with 
H ~ H', we can always successively raise or lower rows of H in a way that 
decreases \th — th>\- 

The first observation is that if H £ S^, then we can choose any square / in G 
on which th is too high (too low) such that / is higher (lower) than three or four 
of its neighbors. Then we lower (raise) the maximal (minimal) row beginning 
with /, thereby decreasing \th — t~h' \ by the length of the row. 

Second, suppose H e S^. Then because H' does not contain any boundary 
vertices, there exists at least one element of B{H) on which t h and th' do 
not agree; assume without loss of generality that there is an / e B(H) with 
T fi(f) > TH'(f)- Then there must be at least one square / in B(H) that is 
maximal along the boundary. If / is higher than three or four of its neighbors, 
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we lower the maximal row beginning with /; if the length of the row is odd, the 
result is in S^. If the length is even, the result remains in . 

We have to be a little more careful in the event that / is higher than two 
neighbors which are boundary faces and lower than the other two, as in the 
following diagram: 



In this case, there is a maximal row beginning with f\ (the face above / 
in the diagram) extending upward and another beginning with f_i extending 
downward, both of which must be lowered before we can decrease the value of 
th on / itself. 

Now, if fi or /_i is the beginning of an even length maximal row, we can 
pull that down without changing the number of cycles in H. If either square is 
the beginning of an odd length maximal row of length greater than one, we can 
first pull up one of the squares adjacent to / on the left or right in the diagram 
(which joins cycles and moves H into S^) and subsequently lower the maximal 
row beginning with f± or /_i (which moves H back into S^). Raising the square 
on the left or right of / is a step in the wrong direction — increasing \th — t~h' 
by one — but it allows us to take three or more in the right direction. 

Suppose, however, that both /i and /_i are local maxima. Then raising 
a square to the left or right of / and subsequently lowering one of /i or /_i 
does not change the value of \th — th'\- However, suppose the longer of the 
two cycles contains j vertices and the other contains k; then we can first raise 
one of the squares to the left or right of / and then lower whichever of f\ or 
/_i belongs to the smaller of the cycles; this has the effect of decreasing the 
length of the smaller cycle. Since the length of the smaller cycle cannot decrease 
indefinitely, if we repeat this process, we must eventually reach a position at 
which it is possible to decrease the value of \th — th'\- □ 

Similar arguments can be made for the connectedness of (Sfj, U S^ +1 ) (~l "K 
under these operations. 

Now that we know that (S^ U S^) fl "K is a connected graph, we can take a 
random walk on this graph. To do this, we define a transition step as follows: 

Choose a face / at random. If H G S^, then if / is higher (lower) than 
three or four neighboring squares and it is possible to lower (raise) the maximal 
(minimal) row beginning with /, then do so. Otherwise do nothing. 
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Figure 13: 
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If H € 5?, then if / is in B{H) and / is higher (lower) than three or four 
neighboring squares and it is possible to lower (raise) the maximal (minimal) 
row beginning with /, then do so. Otherwise, do nothing. If / ^ B(H), but / 
is the beginning of an even length maximal (minimal) row such that lowering 
(raising) that row leaves H in S^, then do so. Otherwise, do nothing. 

Because this transition process describes a simple random walk on a graph 
in which each H has (counting self- loops) degree it has uniform stationary 
distribution. Hence, applying this transition process enough times allows us 
to sample randomly from the space 5j U Sj. The drawback is that we do not 
yet have rigorous bounds on the rate of convergence of this Markov Chain. 
If fl !K is a reasonable fraction of (S^ U we can perform repeated 

Markov searches until we find an element of S^; this gives us a way of choosing 
a Hamiltonian cycle of G at random. A second drawback is that it is not known 
how "reasonable" a fraction fl "K occupies of the entire space (S^ U S^) fl IK. 

Nonetheless, we suspect that the fraction and the convergence rate are in- 
verse polynomial and polynomial respectively in the number of faces of G; we 
have shown this to be the case for some special families of grid graphs. 

7 Almost Taut Examples 

A chipped Aztec diamond of size 2n is created as follows. Start with an Aztec 
Diamond graph of radius 2n. Then successively remove "chips" from the upper 
half of the diamond; here, a chip is a horizontally adjacent pair of vertices 
vi , V2 with a white vertex on the left such that both vertices are incident to 
the unbounded face and none of the vertices above v i and t>2 remains in the 
graph. One can prove that an Aztec diamond of size 2n with k < n chips can be 
partitioned into at fewest n — k disjoint cycles. (The proof follows the lines used 
in the proof of the correctness of our algorithm: that is, we show inductively 
that we can separate the vertices into n distinct cycles with k fixed pairs in a 
configuration satisfying the three properties used in our proof.) Hence k must 
be at least n — 1 in order for the diamond to have a Hamiltonian cycle. We will 
assume in the remainder of this discussion that k is exactly n — 1 . The following 
is an example of a Hamiltonian cycle on a chipped Aztec diamond with n = 4. 

The values of th on the lower half of the diamond are all fixed. Now, label 
the squares with coordinates so that the center squares is (0, 0) and consider 
the 4n — 1 vertical columns of squares in the upper half, each beginning with 
a square of y-coordinate zero. In a sequence of squares jo,-- - , fi from the 
bottom of such a column to the top, the number k for which /j. -< fk-i is 
determined by the boundary data. Let a#(z, j) be the value k for which such 
an event occurs in the ith column for the jth time. The set of i and j for 
which au{i,j) is defined is independent of H. In the example, we have colored 
in the three squares in the upper half on which r is higher than it is on the 
square immediately to the north. Note that the middle square cannot take a 
value higher than than the outer two. That is, a#(0, 1) is less than or equal 
to o/f(l,l) and a#(l,0). We can say that the set of possible H corresponds 
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Figure 14: Chipped Aztec Diamond 



to the set of labelings of an L-shaped Young tableaux with three squares with 
numbers between one and five such that the numbers are non-decreasing in each 
row or column. The reader may deduce a more general correspondence between 
chipped Aztec diamond cycle partitions and Young tableaux labelings. There 
is also a more subtle correspondence between Hamiltonian cycles of a chipped 
Aztec diamond and standard labelings of Young tableaux. 
This next example is even simpler than the last. 

Here there are only three squares on which th takes a higher value than 
it does on the square immediately to the north, and each necessarily occupies 
on of the three columns the squares are in. The reader may show that the 
set of Hamiltonian cycles corresponds to the permutations of three elements 
and deduce a more general result. It is a worthwhile exercise to work out the 
relative sizes of and and try to put upper bounds on the Markov chain 
convergence rates for these simple examples. Intuitively, one would expect that 
in examples that are not so almost-taut, a randomly chosen element of 
would more likely to have one and perhaps several boundary squares beginning 
maximal or minimal rows of odd length (unlike the above two examples, in 
which most would have even length). Hence, one might expect the difference 
in sizes between and to be less pronounced in these cases. It is not yet 
known how to make this rigorous. 

8 Conclusion and Open Questions 

The use of height functions for Hamiltonicity computation and sampling prob- 
lems is a new field of research; there are doubtlessly many improvements to 
be made, empirical studies to be done, bounds to be determined, and general- 
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Figure 15: Low Slack Tower 
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izations to be drawn. There also remain many open problems regarding mean 
properties of restricted difference height functions and performance of Coupling 
from the Past algorithms. We list just a few possible directions for further study. 

1. Can we bound the Markov chain convergence for the set of Hamiltonian 
cycles of a grid graph? Can we generalize the Hamiltonian cycle counting 
algorithms and bounds given in g, p7] , |2S| ], and |3(| for rectangular 
grids? 

2. Can the log of the number of Hamiltonian cycles in a large simply con- 
nected grid graph be approximated with an entropy integral the way the 
number of restricted vertex degree subgraphs can be (see 

3. A close look at the grid graph algorithm shows that we have actually 
produced a characterization of grid graphs that have cycle covers and fail 
to be Hamiltonian; they are those for which "fixing" some r edges in 
H pulls the height function so taut that all the cycle boundary squares 
separating H into a set of m + r regions are held in place (for some m > 0). 
This suggests there may be a faster algorithm for recognizing, for example, 
polyomino graphs with this property. Is there a linear or 0(|.F|Zo<7|.F|) 
algorithm for recognizing Hamiltonian polyomino graphs? 

4. If we assign a weight to each edge in G, there are well-known constrained 
network optimization algorithms that will enable us to find an H parti- 
tioning G into cycles such that the sum of the edges in H has minimal 
weight. What can we say about the problem of finding a Hamiltonian 
cycle of minimal weight? Will the techniques of this paper help in this 
special case of the traveling salesman problem? 

5. Is there a fast algorithm for computing the maximal degree of elements of 

(In some sense, this maximal degree quantifies the notion of "slack" 
in the height function.) 

6. Can we use techniques similar to those for grid graphs to determine the 
Hamiltonicity of hexagonal graphs? (It seems at least plausible; there 
are nice analogs to maximal and minimal square rows. The tricky issue 
is parity — a ^-transformation always changes the number of cycles by 
zero or two. Thus, to decrease the number of cycles, one must find a 
^-transformation that will join three cycles together into one.) If I can 
compute the Hamiltonicity of hexagonal graphs with these techniques, can 
I go still further? Is there a larger class of graphs for which these or similar 
techniques give polynomial algorithms for determining Hamiltonicity or 
other graph invariants? 
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